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ABSTRACT 

We analyze the transit timing variations (TTV) obtained by the Kepler mission for 22 sub-jovian 
planet pairs (17 published, 5 new) that lie close to mean motion resonances. We find that the TTV 
phases for most of these pairs lie close to zero, consistent with an eccentricity distribution that has a 
very low root-mean-squared value of e ~ 0.01; but about a quarter of the pairs possess much higher 
eccentricities, up to e ~ 0.1 — 0.4. For the low-eccentricity pairs, we are able to statistically remove 
the effect of eccentricity to obtain planet masses from TTV data. These masses, together with those 
measured by radial velocity, yield a best fit relation M ~ 3M®(ii/ J?®). This corresponds to a constant 
surface escape velocity of ~ 20 km/ s. 

We separate the planets into two distinct groups, "mid-sized" (those greater than 3i? ffi ), and "com- 
pact" (those smaller). All mid-sized planets are found to be less dense than water and therefore 
contain extensive H/He envelopes, likely comparable in mass to that of their cores. We argue that 
these planets have been significantly sculpted by photoevaporation. Surprisingly, mid-sized planets, a 
minority among Kepler candidates, are discovered exclusively around stars more massive than O.8M0. 
The compact planets, on the other hand, are often denser than water. Combining our density mea- 
surements with those from radial velocity studies, we find that hotter compact planets tend to be 
denser, with the hottest ones reaching rock density. Moreover, hotter planets tend to be smaller in 
sizes. These results can be explained if the compact planets are made of rocky cores overlaid with 
a small amount of hydrogen, < 1% in mass, with water contributing little to their masses or sizes. 
Photoevaporation has exposed bare rocky cores in cases of the hottest planets. Our conclusion that 
these planets are not water-worlds contrasts with some previous studies. 

While mid-sized planets most likely accreted their hydrogen envelope from the proto-planetary disks, 
compact planets could have obtained theirs via either accretion or outgassing. The presence of the 
two distinct classes suggests that 3i? e could be identified as the dividing line between 'hot Neptunes' 
'super-Earths.' 

Subject headings: KOI 137, KOI 157, KOI 168,KOI 244, KOI 870, KOI 952, KOI 1102, KOI 148, KOI 
152, KOI 156, KOI 248, KOI 829, KOI 1270, KOI 1336, KOI 500, KOI 775, KOI 
841, KOI 898, KOI 1215, KOI 1241, KOI 1589, Kepler 11, Kepler 18, Kepler 23, 
Kepler 24, Kepler 25, Kepler 28, Kepler 32 



1. INTRODUCTION 



The spectacular success of the Kepler m ission ( Borucki 



et al. 2010a 2011 



2012 hereafter B12) 



Batalha et al. 

opens our eyes to the world of l ow-mass planets (a lso see 
radial velocity discoveries, e.g. Mayor et al.||2011 |. The 
Kepler mission has uncovered a surprising abundance of 
such planets close to their host stars and to each other. 
While some planets are inferred to have high densities 
and therefore are terrestrial-like, some have low densities 
indicating the presence of a substantial gaseous envelope 
(see references in Table [3| . How do these low-mass plan- 
ets form and migrate? vVhat gives rise to the range in 
planet radius and density? What is the internal com- 
position of these planets? More measurements of planet 
density may help resolving these puzzles. In particular, 
ascertaining whether Kepler planets contain a significant 
amount of water would shed light on the site of their 
construction (see Lopez et al.|2012 for a detailed discus- 
sion) . 

In addition to masses, the planets' eccentricities and 
inclinations are valuable clues. Relative inclinations in 
Kepler systems are inferred to be only a few degrees (e.g., 



Tremaine fc Do ng 2012] |Fabrycky et~aT1[20T2l |Figueira| 
et al.||2012p . This flat configuration invokes images of 
prot o-planet ary disks or planetesimal disks. Eccentrici- 
ties, on t he other hand, are rela tively poorly constrained 
(see, e.g. Moorhead et al.pOll ). Since eccentricities can 
be excited by planet interactions and damped by planet- 
disk interactions or tidal interactions with the host stars, 
they are a fos s il rec ord of past dynamical events. In 
Lithwick et al. ( 2012[ ) (hereafter Paperl) we analyze the 
TTV (transit time variation) signals of six planet pairs 
and show that the TTV phases are consistent with most 
planets having in general small eccentricities. More ec- 
centricity determinations will help us rewind the clock 
and infer the dynamical past of these planets. 

With these two goals in mind, we analyze here the 
TTV signals for pairs of transiting planetary candidates 
that lie near first-order mean-motion resonances (M MR) 



using the Q0-Q6 Kepler table of transit times in Ford 
et al.| p"012| . We also incorporate some results from |Xie' 
2012p . We then ap ply our analytical TTV expressions 
derived in Paper l[ ) to infer planet properties. We in- 



troduce a new element in this work by showing how true 
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planet density can be inferred from the nominal density 
measured by TTV, by using statistical information from 
a large sample. We outline our sa mple selection crite- 
ria in q2] and discuss our r esults in ppl & H But first 
we recap some results from |Paper I| that are of relevance 
here. 

For a pair of transiting planets that lie near a first 
order mean-motion resonance, mutual gravitational per- 
turbations cause each planet's transit timing to vary si- 
nusoidally with a sizable amplitude. The periods, ampli- 
tudes, and phases of the sinusoids encode physical prop- 
erties of the planet pair: 

• The periods of the two sinusoids are the same and 
we call them the super-period: 



= 



\j/P' 



1 
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where P and P' are orbital periods of the inner 
and outer planets, respectively, j stands for the 
j : j — 1 resonance and A = P'/P(j — l)/j — 1 is 
the fractional distance to resonance. Most Kepler 
pairs with clear TTV signals have |A| of order a 
few percent. 

• The amplitudes ( | | and |V|') are functions of both 
planet mass and orbital eccentricity, roughly as 
(disregarding coefficients of order unity) 
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Here, /i is the mass ratio between the inner planet 
and the star and // that of the outer planet. Zf lee 
is a weighted sum of the complex free eccentricity 
in both planets, 



fZfre 



9 z fr, 



(3) 



where / and g are sums of Laplace coeffici ents with 
order-unity values, as listed in Tabic 3 of Paper I. 
The free eccentricity, Zf lee , is the excess in eccen- 
tricity over what is forced by the resonance ('forced 
eccentricity'). We shall often shorten 'free eccen- 
tricity' to 'eccentricity', since the forced parts are 
mostly small for the pairs we consider (~ 10 -3 ). 
These expressions show that amplitude ratio of 
the two sinusoids is roughly related to the ratio 
of planet masses. They also highlight an inherent 
problem in TTV analysis: mass and eccentricity 
are difficult to disentangle. However, if they could 
be disentangled, TTV could be used to probe or- 
bital eccentricity with a sensitivity of a couple per- 
cent (~ |A|). This is superior to other techniques 
like radial velocity or transit shape. 

We define nominal masses for the two planets based 
on the observed TTV amplitudes, 



^nominal — 
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The true mass for the inner planet is m — 
TOnominai/|l - Z bee / (2g A/3) | , and analogously for 
the outer planet. Hence the nominal masses are 
typically upper limits: they are the true masses 
at zero eccentricity, but usually lie a bove the true 
masses when eccentricity is non-zero ( Paper 1 1 . 



• Phases of the TTV sinusoids encode information 
about the free eccentricity. TTV phases are de- 
fined relative to the time when the longitude of 
conjunction points at the observer. The two phases 
(0ttvj0ttv) should lie at (0, 7r) for zero eccentric- 
ity^] while for most other eccentricities, the two 
phases can take other values but remain largely 
anti-correlated (4>' ttv ~ <fittv + tt)- Equivalently, if 
the phases are close to (0, 7r), then the pair likely 
have small eccentricity. We show below that TTV 
phase can be used to infer planet eccentricity (<|3]), 
and consequently, genuine planet masses (§41). 

2. A SAMPLE OF 22 PAIRS 

Here we describe how our sample of 22 TTV pairs are 
assembled. 

Transit tim es for quarters 0-6 are published by [~ 



et al. (2012) for a large number of Kepler candidates. 
We select from this list all pairs within |A| < 8% from 
a first-order MMR (2:1, 3:2, 4:3, 5:4). For each pair, we 
obtain average orbital periods and TTV amplit udes and 
phases using the fitting procedure described in Paper I| 
We then calculate the periodogram for each TTV series 
and include a pair for further consideration only when 
both planets have periodograms that show a clear ex- 
cess of powen by visual inspection, at the desired super- 
period (eq. fl]). We only include pairs with super-period 
shorter than 1000 days, to ensure that more than half of 
the TTV sinusoids have been observed. Amplitudes and 
phases of the TTV sinusoids are measured at this super- 
period, using a least-squares fit (Paper I). Weeding out 



pairs whose sinusoid phases are not roughly out of phase 
with each other, we find a total of 20 pairs that pass 
these thresholds, 6 of which are the confirmed systems 
analyzed in Paper I. 

Compared to the exquisite TTV d ata for koi 137.01/02 
(Kepler-18c/d, [Cochran et al.|2011[ ), TTV data for typi- 



cal Kepler candidates are less accurate, since planets are 
typically smaller in size and sampled at a lower cadence. 
Our stringent weeding criteria are necessary to overcome 
random noise, and to isolate perturbations from other 
planets. The latter is a prominent issue in Kepler data, 
as planet pairs are often accompanied by other planets 
with near resonant period ratios (e.g., KOI-82, KOI-500, 
KOI-2169). 

In Fig. |11[ we present the sinusoidal fits, error- 
bars, and periodograms for 6 of the above pairs, koi 
156.01/03, koi 775.02/01, koi 841.01/02, koi 1215.01/02, 
koi 1241.02/01. Exce pt for koi 157.01/02 (Kepler-11 e/f, 
Lissauer et al. 2011a ), these pairs have not been identified 
in the liter ature before. For pairs previously identified by 
Xie| ( |2012[) , we adopt his TTV fits for these systems as 
they are obtained using a longer baseline (through quar- 
ter 9). We also include 2 new pairs from that work, koi 
829.01/03, and koi 1336.01/02, that our shorter dura- 

1 Prom now on we call the inner planet's 0ttv the TTV phase. 
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Fig. 1. — Nominal planet density versus radius for the TTV 
sample in Table [TJ Nominal densities are typically higher than 
true densities. Magenta stars stand for previously con firmed pairs 
and green stars those identified in this study (Fig. \U\ . Error bars 
mark 68% confidence limits on the nominal densities (solid error 
bars for inner planets and dashed for outer). Typical uncertain- 
tics in radius are ~ 10% and are not displayed here. The TTV 
method is currently sensitive to planet masses as small as 2M($. 
The colored curves, marked with 'Iron', 'Rock' and 'Water' respec- 
tively, represe nt the theoretical ma ss-radius relation for these pure 
compositions | |Fortney et al,|2007| ). 

tion data failed to pick up. But we do not use 3 objects 
from that paper, koi 869.03/02, koi 877.01/02 and koi 
880.01/02: the first two because their TTV phases are 
not roughly anti-correlated, which may indicate pollu- 
tion, the last one because of an unphysically large den- 
sity (and mass) for the inner planet: it requires a Jupiter 
mass for a 2.7i?© planet. It is also likely contaminated 

by a nearby planet. 

So together with the 6 planet pairs analyzed in Paper 
I ( we have a total of 22 convincing TTV pairs (5 new, 17 
published). TTV solutions for these pairs, including their 
nominal densities (nominal mass divided by volume), are 
listed in Table [T] and depicted in Fig. [T] 

2.1. Selection Effects 

There are 151 planet pairs in the |B12| catalog that lie 
within |A| < 5% from a first-order MMR. Our sample 
of 22 pairs is ~ 15% of the total. Since the strengths 
of TTV signals depend on planet mass and eccentricity, 
it is of concern whether our sample is biased towards 
planets that have higher density or higher eccentricity 
than average. Moreover, we sample only planets near 
MMRs, which may in principle have different properties 
than those far from resonance. 

We argue here that, in comparison with the B12 cata- 
log, our sample favours larger planets, as well as planets 
that have suitable orbital periods, but not necessarily 
planets with higher density or eccentricity. First, larger 
planets show deeper transits and allow more precise tim- 
ing. The median size of the B12 pairs is 2.18i?m while 



bital periods and hence long super-periods. It also ex- 
cludes planets too close to their stars, which have smaller 
intrinsic TTV amplitudes and less accurate TTV mea- 
surements due to their shorter transit duration. The me- 
dian inner planet period in our sample is 8.04 days, with 
a standard deviation from the median of 5.98 days. In 
contrast, values for the |B12| sample are 7.88 days and 
14.1 days. This means that we are mostly probing pairs 
clustered around ~ 8 day orbital periods. 

These two biases alone could explain why our sample 
is only 15% of the overall sample of candidate pairs. In 
fact, when restricting to pairs similar in properties to 
our sample (inner periods between 5 and 20 days, planet 
radii g reate r than 1.7i?®, |A| < 5%), we find 48 pairs 
in the B12 catalog. This leaves little room to suspect 
that we are significantly biased toward higher density 
and eccentricity. 

A related issue concerns sensitivity in nominal mass. 
Currently Kepler can probe TTV signals with amplitude 
upward of ~ 5 minutes. Adopting |A| ps 2%, and an in- 
ner period of 8 days, this implies that TTV's are sensitive 
to planets with nominal masses upward of 5-M© , around a 
sun-like star. Sensitivity around lower-mass stars is even 
better, ~ 2M ffi (Fig. nl), thanks to their deeper transit 
at the same planetary size. Lastly, we find that, with a 
factor of 2 longer observation, or moderate improvement 
in transit timing precisions, one should be able to detect 
TTV sinusoids in many more pairs than reported here. 

We argue that our near -reson ant pairs are not special. 
First, planet pairs in the |B12| sample show similar size 
distribution, whether they are close to an MMR or not. 
Moreover, while our TTV sample is restricted to planets 
near first-order MMRs, the RV planets we include in our 
analysis below are not, and those appear to have similar 
properties. 

3. MEASURING PLANET ECCENTRICITIES 



One of the most important results from Paper I is the 
utility of the TTV phase (0ttv) for inferring the value of 
the free eccentricity. If a significant fraction of systems 
have 0ttv ~ 0, it can be inferred that most of those have 



Z free | < |A|. In Paper I we find that the TTV phases of 



our sample has a significantly larger median of 2.72i? ffi . 
Second, our TTV sample excludes planets with long or- 



the six analyzed pairs are clustered near zero, indicating 
that in general the free eccentricity is small (|Zf ree | < A). 

Now armed with 22 pairs, we explore the eccentricity 
distribution further. We find that for this larger sam- 
ple too the TTV phases are not uniformly distributed, 
but cluster around zero (Fig. [2]), again indicating that 
the general population possesses small free eccentricity. 
However, there are a number of pairs that buck this 
trend. These are marked in Fig. [3j based on three in- 
dices: TTV phase, nominal density, and transit duration. 
Both the TTV phases and nominal densities (compared 
to other planets of similar sizes) for Koi 775, 1241, 1270 
are large, suggesting high eccentricities. Since for the 
high-e population, cj> ttv is uniformly distributed between 
(— 7T, 7r), we expect to find pairs that have high-e but low 
TTV phases. And indeed, pairs Koi 152, 829, 1215, while 
showing small phases, have high nominal densities. The 
transit durations for these planets are often longer than 
expected for circular orbits. We categorize these 6 pairs 
as the high-e population and analyze them in |5.1| It is 
interesting to note that all planets with nominal densities 
above ~ 15g/cm 3 are now classified as high-e. 
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TABLE 1 

TTV ANALYSIS OF 22 PLANET PAIRS 
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Note. — A list of our sample of 22 TTV pairs. We present their measured TTV amplitudes and phases, resulting nominal density 
(nominal mass divided by volume, in unit of g/cm 3 ), and other relevant parameters. P and P' are orbital periods in days, R and R' are 
planet radii in unit of earth Radius (R&) The first 8 pairs are confirmed systems (corresponding to Kepler 18c/d, llb/c, lle/f, 23 b/c, 25 
b/c, 28 b/ c, 32 b/c, 24 b/c, respectively), six of which were analyzed in |Paperl| TTV parameters for the next 8 pairs are adopted from 
|Xie| | |2012[ |, and TTV sinuso ids for the last 5 pairs are reported here for the hrst time (Fig. We adopt stellar parameters and planetary 
radii trom the KIC catalog l |Brown et al.|2011| , except for the confirmed systems where more refined parameters may be available. Pairs 
marked with the * sign are suspected to nave significant free eccentricity (jj3|. Typical error bars in planetary radii are of order 10%, 
smaller than those in TTV amplitudes. These arc not included in calculating uncertainties in nominal densities. 



The remaining 16 pairs (~ 75%) are consistent with 
low eccentricity. We assume them to be a homogeneous 
population satisfying a single eccentricity distribution. 
TTV phases depend on Zf veo , a weighted sum of the two 
planets' free eccentricities (Eq. [3]). Splitting into real 
and imaginary parts, 



(5) 



and adopting a 1-D gaussian distribution for these com- 
ponents (for both the inner and outer planets) , P{e x ) = 
Pa{e x ), P{e y ) = P<j{e y ), with 



P*(e x ) = 



1 



2na 



exp(- 



(G) 



of 



TTV 

that 



we can calculate the resulting distribution 
phases. We assume e x and e y are uncorrelated 
the phases of the complex eccentricities are r andom, a s 
they would be due to secular or GR precession ( Paper I ) . 
TTV phases are either uniformly distributed or clustered 
around zero, depending on the relative ratio of cr/|A|, 
and on the resonance of relevance. Generating a mock 
catalogue that has the same period ratios as our low-e 
sample, we produce different phase distributions for dif- 
ferent values of a (Fig. El. The current sample lead us to 
conclude that a value of a « 0.007 provides the best fit, 
with an uncertainty in a of order 20%. This value cor- 



responds to a RMS value in eccentricity (e = , - , . „ 

of 0.01. So it appears that the low-e population has free 
eccentricities that are of order (but somewhat smaller 



than) the typical resonance offset (|A| ~ 0.02). 

This exercise demonstrates the exceptional sensitivity 
of TTV measurements to even small values of eccentric- 
ity. We discuss the implication of these results in ^5.1 



4. TRUE DENSITIES 

Employing our newly found eccentricity distribution, 
we perform statistical simulations to correct for the ec- 
centricity effect, and thereby obtain true planet densities 
for the low-e population, subject to statistical and mea- 
surement errors. Densities of the high-e sample are not 
recoverable. 

For each planet pair, we generate complex free eccen- 



tricities (Zf, 



) weighted by the intrinsic distribution 



ree > *-f ree 

inferred above, with a — 0.007 (Eq. ph. Each pair of ec- 
centricities produces correction factors relating the true 
and nominal masses (see below Eq. El). Keeping only 
those generated eccentricity-pairs that reproduce the ob- 
served TTV phases (</W, 4>'ttv) within the 68% error bars, 
the median correction factor then gives the median true 
density. This procedure introduces an uncertainty which 
we quantify by the 68% bounds in the distribution of 
correction factors. For the total uncertainty in density 
measurement, we add these bounds quadratically to the 
error-bars in TTV amplitudes. Results are listed in Ta- 
ble [2l The ratio of median to nominal density depends 
on the observed TTV phases, but is typically of order 

2 We have also experimented with other forms of eccentricity 
distribution, e.g., a <5-function plus a Gaussian, 1-D Gaussian, etc., 
but find no improved fit to the data. 
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Fig. 2. — Phases of the TTV sinusoids for the inner (0ttv) and 
outer (<j>' ttv ) planets for our 22 pairs. Pairs with very small eccen- 
tricity (-C |A|) would lie near (0,180°). For higher eccentricity, 
phases should lie along the various grey curves (from 0.25|A| to 
2|A| in increments of 0.25|A|, dark grey for 2:1 MMR and lighter 
grey for 3:2), becoming increasingly evenly distributed along the 
diagonal line, which corresponds to anti-correlated phases. His- 
tograms for the measured </>ttv and (p' ttv are displayed at the top and 
side panels, respectively. The measured TTV phases cluster around 
(0,180°). In this and all subsequent figures, red denotes planets 
around stars more massive than 0.8Mq, and blue those around less 
massive stars. Error-bars indicate 68% confidence limit. Phases of 
three pairs (labelled with their koi numbers) deviate significantly 
from zero. These also exhibit anomalously high nominal densities 
(Table [Tj, suggesting eccentricity > |A|. The pair koi 157.06/01 
(Keplerllb/c) also has high phases (the red point at (97°, 237°)), 
but this likely results from its dominant resonance (5:4) being close 
to other nearby resonances (4:3, 6:5, etc). 

0.5. Even for planets with TTV phases consistent with 
zero, the correction factor falls below unity. This is be- 
cause the most likely eccentricity value lies not at zero, 
but at e ~ a. 

The median densities and their error-bars are plotted 
in Fig. 5] Combining our results with a sample of low- 
mass planets that have previously determined masses 
(Table |3j mostly through radial velocity), we find a best- 
fit mass-radius relation of 



M ps 3M a 



R 
Rm 



(7) 



or a density p ps 3p©(i?/i?©) -2 where the density of the 
Earth is /?© = 5.5g/cm 3 . This best-fit d iffers from the 
one th at applies for solar system planets dLissauer et al. 
2011b[ ), M = {R/R <s ) 2m M ( 
stead 



or p ps p e (R/R m ) \ In- 
it corresponds to one that has a constant surface 



escape velocity of 20 km/ s. This is fortuitously close to 
the sound speed of hydrogen plasma at 10 4 K (13 km/ s), 
suggesting that the process of photoevaporation of hy- 
drogen may be involved in some way. 

In the next section, we discuss further details of our 
findings, as well as their implications. 

5. DISCUSSIONS 
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Fig. 3. — Selecting pairs with high eccentricity by their transit 
duration, nominal densities, and TTV phases. We plot TTV phase 
(4>ttv, top panel) and transit duration (lower panel) versus nominal 
density for all 22 pairs. The transit duration is normalized by 
^dur 0i the duration of a circular orbit with zero impact parameter. 
All circular orbits should lie below the line of unity, while eccentric 
orbits can lie above or below. Three pairs (Koi 775, 1241, 1270, 
marked out in Fig. pi exhibit high nominal densities, large TTV 
phases, and longer tnan unity transit duration. These systems 
almost certainly have high free eccentricities (> a few X A). Three 
pairs (Koi 152, 829, 1215) have low TTV phases. But their high 
nominal densities, and in many cases, longer than unity transit 
duration, suggest that they possess high eccentricity. Together, 
these 6 pairs constitute our high-e population. Transits of koi 248 
last longer than expected of circular orbits. We suspect radius 
determination for the host star is too small by 10 — 30% (but we 
neglect this in the following analysis). 

We first divide planets into two groups: 'mid-sized' 
(those with R > 3i?©) and 'compact' (those smaller). 
This division line is partly inspired by an observation of 
the overall Kepler sample (B12 ): there is a drastic fall-off 
in planet numbers around 3-ft© (right panel of Fig. [7]), 
suggesting that a transition occurs. We aim to discover 
the nature of this transition. 

5.1. Dynamical Excitation 

For the low-eccentricity population, we obtain an RMS 
eccentricity of e ps 0.01. Although small, this value is un- 
expected. To explain the asymmetric pile-up of planet 
pairs on the far side of mean-motion resonances , Lith- 
wick fc Wu| ( |2012[ ) and |Batygin fc Morbidelh| ( |2012| argue 
that Kepler planets have experienced protracted orbital 
damping. This could arise, for instance, during their 
tenure living in disks. If this explanation is correct, the 
small free eccentricities observed in TTV pairs would re- 
quire one to invoke a new source of eccentricity in the 
system to re-excite the planets. 

Eccentricities for the high-e population can be esti- 
mated assuming that the true densities for these pairs 
follow the fiducial mass-radius relation (Eq. 17 1. This 
yields e ps \z frcc \ ps (p n0 m/pfiduciai) |A| ps 0.1 - OA. 

Such high values are puzzling. To understand their 
origin, we look for correlations between eccentricity and 
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TABLE 2 

Density after Correcting for Eccentricity Effect 
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Fig. 4. — The distribution of TTV phases. The black histogram 
is the observed distribution for the 16 low-e pairs. The colored 
curves are simulation results at different values of eccentricity dis- 
persion (Eq. [61, calculated for a mock catalog that has the same 
period ratios as the observed one. The curve with a = 0.007 best 
reproduces the observation, with an uncertainty of ~ 20%. The 
RMS eccentricity of individual planets is \/2cr S3 0.01. 
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Fig. 5. — Similar to Fig. ^ but showing the corrected densities 
for the 16 low-e pairs, after adopting a Gaussian model for the 
eccentricity distribution with a = 0.007 (Fig. H. We apply no 
correction for Kepler llb/c (koi 157.01/02). Also shown here are 
density measurements from Table [3] (golden squares, mostly by 
RV, except for Kepler 36b/c). The RV sample exhibits similar 
behavior to the TTV sample: larger planets are less dense, and 
are comparable in mass to the smaller ones. The points are best 
fit by the thick black line: M = 3M^(R/ Rq), corresponding to 
a constant surface escape velocity of v eBC = 20km/s. All planets 
satisfy d csc > 13 km/ s, the sound speed of 10 4 K hydrogen plasma. 
All densities fall below ~ 15g/cm 3 . 
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Note. — Planet densities after applying the statistical 
model for the eccentricity distribution. The column p lists 
the median density for the inner planet (and primed quanti- 
ties for the outer planet), a product of the median correction 
factor (p/pnom) and the nominal density. Error-bars on the 
correction factor reflect the width within which 68% of all 
possible solutions fall. Error-bars on the median density are 
a quadratic sum of uncertainties in the correction factor and 
uncertainties in nominal mass (measurement error in TTV 
amplitude). The last column lists the average of the two me- 
dian eccentricities, as a rough indicator for the magnitude of 
eccentricity in the pair. Our values for Koi 137.01/02 (Kepler 
18c/d) agree wit h previous determinat ions of p = 0.59 ± 0.07, 
p' = 0.27 ± 0.03 ( |Cochran et al.|2011[ l. We do not perform a 
correction for Koi 157.06/01 (Kepler llb/c) as its TTV phase 
is polluted by other resonances. 



other properties. We do not observe significant correla- 
tion with stellar mass, planet size or planet period. In 
particular, mid-sized planets and compact planets alike 
contain high-e pairs. We do note that all high-e pairs 
reside near the 2:1 resonance. Had they been near 3:2 

TABLE 3 

Literature Determinations of Planet Density 



Planet 



gCge) 



GJ1214b 
GJ3470b 
GJ436b 
Kepler-4b 
Kcpler-lOb 
HAT-P-26b 
HAT-P-llb 
CoRoT-7b 
55 Cnc e 



2.68 
4.1 
4.0 
3.99 
1.42 
5.8 
4.9 
1.6 
2.2 



A/(Mffi) p(g/cm J ) 



6.55 

14.5 

24.3 

24.5 

4.56 

18 

26 

7.4 

8.6 



1.87 ± 0.4 
1.16 ±0.17 
2.09 ±0.14 
1.9 ±0.4 

8.8 ±2.5 
0.4 ±0.1 
1.2 ± 0.13 
10.4 ± 1.8 

5.9 ± 1.3 



kepler-20b 1.9 8.7 7.1 ± 1.8 

kepler-20c 3.1 16.1 3.0 ± 0.8 

Kepler-36b 1.5 4.5 7.5 ± 0.7 

Kepler-36c 3.7 | 8.1 0.89 ± 0.06 

Note. — Transiting planets with mass determinations using 
the radial velocity technique (except for Kepler-36b/c which is 
through TTV). Here, we only include planets that are smaller 
than 9Rm in size. 



Reference 

Charbonncau ct al. (2009) 
Bonfils ct al. (2012) 
Gillon et aCf||2007[) 
Borucki et al| | [2010b[ | 
Batalha ct al. (2011) 
Hartman c f~al.| (|201l) 
Bako s et aff^loT "^ 
Hatzes et af]||201l| 
Endl et al.||20T2|~ 
Winn et ahjl|201l| 
Gautier et at] ( |2012[ 
Gautier etlif| |2012[ 
Carter et aft(|2012t 



Carter ct al. (2012) 
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Fig. 6. — Hydrogen mass fractions on planets today (grey points) 
and those that are lost by 5 Gyr of photoevaporation (colored) , for 
objects in Fig. [5] The grey points assume the planet is composed 
of a rocky core covered by a hydrogen atmosphere of solar compo- 
sition, and the colored points assumes the hydrogen is eroded by 
energy-limited escape (Eq. [81. Given uncertainties in radius mea- 
surement, we are only sensitive to hydrogen mass fraction > 10~' 5 . 
All grey symbols falling below this value should be regarded as 
upper limits. For a given planet, if its grey symbol sits above its 
coloured one, it could have retained much of its hydrogen envelope 
over its lifetime. While if its grey symbol sits lower, as is the case 
for most compact objects, the planet should be a bare core. 

or even closer resonances, they might be unstable. We 
speculate that, while 25% of our pairs have high-e, the 
process that excites the eccentricities to such high values 
could have happened more often, but only those near 2:1 
remain. 

We consider the hypothesis that these high-e popula- 
tion arise from planet-planet scattering. If so, one ex- 
pects that the relative inclinations between planets to be 
high, since close-encounters tend to excite inclinations as 
well as eccentricities. However, Koi 152, 775, 829 each 
show 3 transiting planets^] Such a high incidence of tran- 
siting planets argues against high relative inclination in 
the system, and hence against planet scattering. 

The presence of free eccentricity allows one to place 
a lower bound on the tidal quality factor (Q, e.g. 



Gol- 



dreich k. Soter||1966[ ). Considering the high-e pair koi 
1270.01/02, for example, in order for its inner planet to 
retain eccentricity, it needs to have a weak tidal dissipa- 
tion with Q > 1000, approaching the Q values estimated 
for gas-rich planets in our Solar system (Q ~ 10 5 ). In 
fact, if tidal dissipation in these planets is as efficient as 
that in Earth (tidal quality factor Q w 10), tidal inter- 
action between planets and the host stars should have 
reduced most planets' eccentricities to zero. 

5.2. Mid-sized planets, R > 3i?® 
5.2.1. Structure, Photoevaporation 

3 Intriguingly, the third planets in these systems all lie close to 
2:1 MMR with one of the TTV planets. 



Mid-sized planets are a minority. They comprise 23% 
of the total Kepler planet candidates)^] Figure [H] shows 
that mid-sized planets are less dense than water and must 
have extensive H/He en velopes. Based both on the evo- 



lutionary calculations of Fortney et al 



(2007), and on our 



simple model (see below), we conclude that the gas en 
velopes of these pla nets are comparable in mass to their 



mi 

solid cores (al so see |Gillon et al.|200 7 Bakos et al.||2010 
Baraffe et aL|2008[ for KV planets). 



We calculate the fractional masses in their hydrogen 
envelopes that is required to account for their observed 
densities and sizes. We assume that the cores of these 
planets are made of pure roc k, with the co r e radi i and 
masses related by Eq. (7) of Fortney et al. (2007). We 
assume the photospheres of all planets are at a pressure 
of 1 bar, and a temperature that is solely determined by 
irradiation, T = T oq , where T cq = (R*/2a) 1/2 is the 
equilibrium temperature of a black-body at distance a 
around a star with temperature and radius R*. We 
integrate the thermal structure of a solar-composition 
gas inward t owards the core surfa ce, using the equa tion 
of state from Saumon et al. ( |1995[ ), and opacity from | A1-] 
lard et al. ( 200ip . The temperature structure we obtain is 
nearly isothermal near the surface and adiabatic in the 
deep interior. The core mass is adjusted until the ob- 
served total mass is reproduced. The result depends on 
the assumed internal luminosity since the latter affects 
the temperature gradient and therefore the pressure scale 
height at the base of the envelope. For a luminosity of 
lO~ n £0 (the current value on Neptune), the mid-sized 
planets require hydrogen envelopes of order 2 — 30% in 
mass (Fig. [6]). The fractions rise to 10 — 90% when the 
internal luminosity drops by a factor of 100. 

The mid-sized planets exhibit surface escape velocities 
comparable to the sound speed of an HII region (10 4 K, 
13km/s), suggesting that photoevaporation is important 
for their evolution. This is corroborated by the absence 
of mid-si zed planets with periods shorter than a few days 
(see Fig. 10). Here, we investigate this possibility quan- 



titatively. 

We estimate mass-loss due to photoevaporation follow- 
ing Murray-Clay et al. (2009) for XUV photoevapora tion. 
X-ray will enhance the rate (Owen & Jackson 20121 but 
is not considered here. Let the stellar XUV luminosity be 
L\jy . When this flux is absorbed by the planetary atmo- 
sphere, it drives a thermal expansion and outflow. In the 
limit of low flux, the outflow rate can be approximated 
by the so-called energy-limited formula, 



Mi 



L\jvttR p 
Ana 2 v 2 



1.2 x 10 10 g/s 



Luv 
5 x 10" 6 L 



o 



0.1AIA 2 / R p 



\5Ra 



3 nSMg 



(8) 



where we have scaled the XUV luminosity by 5 x 
10~ 6 L Q , roughly the current solar value, and have as- 
sumed an efficiency (e) of unityQ We have verified that 



4 Here, we only include candidates in systems with multiple tran- 
siting objects, as they less polluted by false-positives l |Lissauer et al.| 
|2012| . 

S— We also equate the absorption cross section to the size of the 
planet disc because the atmosphere pressure scale height is very 
small, H/Rp fts 0.01 {R P /2R <S ) (15M & /M P ) (T/1000K) < 1. 
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the radiation/recombination- limited mass-loss rates, rel- 
evant when heat loss is significant (Murray- Clay et al. 
20091, do not apply: they lie above the energy-limited 
rate for all our objects (compact and mid-sized). We 
further no te that the fo rmally defined sonic radius, r s = 

frequently lies inside the physi- 
surface-escape planets we consider 



or the 



|1965| ) 

low su 



GM p /2<? s (|Parker| 
cal radius 
here. 

Estimates for the fractional mass loss, integrated over 5 
Gyrs, are presented in Fig. [6j We have uniformly taken 
the XUV luminosity to be 5 x 10 -6 of the bolometric 
value, and have adopted values for stellar r adius and ef- 
fectiv e temperature from the KIC catalogue (Brown et al. 
20111 for the TTV planets, and from discovery papers 
for planets in Table [3]) j^] We find that mid-sized plan- 
ets could shed between a few percent to a few tenths of 
their masses, comparable to the hydrogen mass fractions 
inferred to be present on these planets (grey symbols in 
Fig. [6]). Our estimates for the mass-loss rate are for the 
current systems. The rates were likely higher in the past 
when the stars were more ch romospherically active and 
brighter in the FUV(see, e.g. |Lammer et al.p 



2004 Ribas 



et al.|200"5 ), and also when the plan ets were younger and 



nuttier (also see Lopez et al. 
tional XUV flux may be hig 
than the Sun 



2012). Moreover, the frac- 



ler for stars lower in mass 
We therefore conclude that photoevapo- 
ration has significantly sculpted these mid-sized planets. 

In the following, we further argue that photoevapora- 
tion may explain the mass-radius relation for mid-sized 
planets (eq. [7|, and naturally produces planets with 
co mparable core /env elope masses (also see discussions 
Lopez et aL|2012| ). 



-motoevaporation proceeds healthily as described by 
Eq. Q as long as the surface escape velocity of the 
planet remains low. When the latter rises, the sonic ra- 
dius (for a wind of ionized hydrogen) moves much beyond 
the planet surface. This throttles the photoevaporation 
rate as the radiation/recombination- limited mass-loss, 
now relevant, decreases exponentially with the expan- 
sion of the sonic surface (Murray-Clay et al. 1120091). Now 



imagine a gas-rich planet that was exposed to the strong 
XUV flux from its young host star before it has ther- 
mally contracted, either because it was for med in-situ or 
was m igrated very quickly (eq. ([8]), also see |Rogers et al. 
2011 1. If such a planet started with much more gas than 
that in its core, evolutionary calculations ( |Fortney et al. 
20071 show that its size will expand during mass loss, 
leading to further evaporation. This process continues 
until mass of the solid core contributes significantly to 
gravity. At this point, the planet shrinks when it loses 
mass until the escape speed from its surface exceeds the 
thermal speed of ionized hydrogen. Planets thus sculpted 
are expected to have envelopes that are lower or compa- 
rable in mass to their cores. These arguments highlight 
the need for a more accurate model, where one accounts 
for the time evolution of both the XUV luminosity and 
planet structure. 

Regardless of the evolution, the mid-sized planets al- 
most certainly acquired their massive hydrogen envelopes 
by accreting from their proto-planctary disks. 



6 Some of the stars have only recently become sub-giants but we 
ignore this complication here 
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Fig. 7. — Size of Kepler planets versus their host mass (left panel) 
and distribution of planet sizes (right panel). Here, we include 
only planet candidates from the B12 catalog that are in multiple 
systems (black and colored points, where the latter are the sample 
considered in this paper). The black histogram is the total size 
distribution, and the colored curves are break-downs into different 
stellar mass bins. Overall, there is a precipitous drop of planet 
numbers at the dividing line at 3-R®. And mid-sized planets are 
absent around stars less massive than 0.8Mq (with the exceptions 
of GJ 436b and GJ 3470b, both transiting planets discovered by 
ground-based telescopes). As the stellar mass increases from 0.4 
to 1.4Mq, the fraction of mid-sized planets rises from to unity. 

5.2.2. Correlation with Stellar Mass 

We focus here on one remarkable correlation between 
mid-sized planets and their host properties, discovered 
in the published Kepler catalogue. In Fig. [7J we display 
the distribution of planet sizes as a function of their host 
mass. Mid-sized planets occur almost exclusively around 
stars more massive than 0.8M Q , or T g > 5000 K on 
the main sequence. Quantitatively, the fraction of mid- 
sized planets, 23% in the overall sample, is nearly zero 
for stars less massive than 0.8M Q , and of order half for 
stars more massive than l.lM©. By contrast, compact 
planets are more equitably distributed around stars of all 
spectral types, albeit with a trend that average planet 
sizes decrease around lower mass stars. 

The fact that K/M-dwarfs do not host mid-sized, gas- 
rich planet^] has not been previously noticed. We offer 
two possible explanations but no final resolution. First, 
photoevaporation around lower mass stars could have 
been more severe than around F/G dwarfs as a result 
of more potent and more prolonged chromosphcric ac- 
tivity. Second, planetary cores may take longer to form 
around lower mass stars, so that they cannot accrete gas 
before the protoplanetary disks disperses. Both scenarios 
can also explain why compact planets tend to be smaller 
around lower-mass stars. 

5.3. Compact planets, R < 3i?® 

We turn now to the compact planets. The mea- 
sured densities of compact planets range from ~ 1 to 
~ 10g/cm 3 (Fig. [5]). The largest compact planets 
(R ~ 2 — 3i?®) have densities that hover around the 
water line (p ~ 2g/cm 3 ). At any given size, planets 
span almost a decade in density. While the mid-sized 



7 |Gaidos et al.| ( |20"l2t argued that some K/M-dwarfs in the KIC 
are misclassmed giants. If true, the sizes of their planets would 
increase to mid-sized. However, this would also mean only sub- 
giants have mid-sized planets, which would be puzzling. 
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Fig. 8. — Hotter compact planets are denser. We scale the den- 
sity of each compact planet to that of a water sphere of the same 
mass (left-hand ordinate), and to a rocky sphere (right-hand or- 
dinate). Although the densities of watery and rocky spheres are 
mass-dependent, the conversion factor between the two is roughly 
constant (Fig. pi. The horizontal axis is the black-body equilib- 
rium temperature of the planet. We include all planets in Fig. |5] 
that satisfy R < 3.2R^. While cool planets are roughly compatible 
with being water-worlds, the hottest ones have densities of rock. In 
contrast to density, planet masses show no systematic dependence 
on equilibrium temperature. 

planets are most likely cores overlaid with massive H/Hc 
envelopes, the structure of compact planets is less cer- 
tain. They can either be water-rock mixtures, or rocky 
cores enveloped by a small amount of hydrogen. It is 
important to resolve this degeneracy because a water- 
laden planet would favour the scenario where the planet 
was constructed outside the ice-line and then migrated 
inwards; a rocky core would be consistent with in-situ 
formation. Although density measurements of individ- 
ual planets cannot break this inherent degeneracy (see, 



e.g., |Adams et al. 2008[ ), we attempt to break it with 
our statistical sample. Fortuitously, the sample straddles 
the region where hydrogen photoevaporation can signif- 
icantly alter planet sizes. 

5.3.1. Density Correlates with Temperature 

Planet densities exhibit a remarkable correlation with 
environment: planets with hotter equilibrium temper- 
atures are in general denser (Fig. M). In that figure, 
we scale each planet's density to the density of a pure- 
water (or pure-rock) planet of the same mass (see the 
theory curves in Fig. [BJj^] Cold planets (T cq sa 600 K) 
are compatible with being pure water, while hot plan- 
ets (T eq > 1500 K) are compatible with being pu r e rock . 
Similar behavior has been noted by |Lopez et al.| ( |2012[ ). 

A second piece of evidence for this correlation is tur- 
nished by studyingthe density ratio within planet pairs 
(top panel of Fig. pi). Unlike individual planet mass, the 
mass ratio of two planets within a pair can be largely de- 

8 If we do not scale the densities, the correlation is still present 
but is more scattered. 
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Fig. 9. — Properties of planet pairs plotted versus the ratio of 
their orbital periods. The top panel shows the ratio of nominal den- 
sity between the inner and outer planets for the compact planets 
in our sample. For pairs that lie close to the 2:1 MMR, the in- 
ner planets are invariably denser, while pairs around more closely- 
spaced MMRs appear indistinguishable from each other. The bot- 
tom panel shows the size r atios for these pairs (colored points), as 
well as for all pairs in the[B12 catalog that satisfy P < 10 days 
and where both components are compact (grey triangles). The in- 
ner planets tend to be smaller, more strikingly so when the planet 
pairs are spaced further apart. 

termined without knowledge of the free eccentricity. So 
here we simply plot the density contrast as the ratio of 
their nominal densities. We find that the inner planets 
in pairs near 2:1 resonance are always denser by a factor 
of a few compared to the outer ones, while pairs near 
closer resonances (5:4, 4:3, 3:2) have more comparable 
densities]^] So some process appears to be at work to 
systematically compactify the inner planets. 

The last piece of evid ence c omes from studying sizes of 
compact planets in the |B12| catalog. The bottom panel 
of Fig. [9] shows the size ratio within a planet pair as a 
function of their period ratio. We only plot pairs with 
inner period P < 10 days. Inner planets tend to be 
smaller than their companions. This trend is more sig- 
nificant for more widely spaced pairs. And it is absent 
for pairs orbiting too far from the star (P > 10 days). 
We also plot planet sizes versus orbital period or equilib- 
rium temperature in Fig. |10| One observes that compact 
planets further away from tneir stars ca n be larger than 
those closer in (also see, e.g., Fig. 8 of B12). The me- 
dian size for planets inward of 5 days is R ~ 1.5i?®, 
compared with R ~ 2.5i?® for planets between 5 and 20 
days. These results suggest that planets inward of ~ 10 
days, or ~ 1000 K around a sun-like star, show signs of 
compactification. 

5.3.2. Rock or Water? 

A naive explanation for Figs. [8|9] is that all compact 
planets started out as water/rock mixtures with water 

9 One notable exception is Kepler-36 where the inner planet 
is 8 times denser than the ou ter one, while the two are near 6:7 
resonance | |Carter et al.|2012[ . 
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Fig. 10. — Planet sizes in the |B12j sample, plotted versus orbital 
period (above) and black-body equilibrium temperature (below). 
We only include systems with multiple candidates. Red points 
indicate those around more massive stars, and blue those around 
K & M dwarfs. 

being the dominant component, but the hottest ones 
have experienced total water removal and became rocky 
spheres. 

This, however, seems difficult on theoretical grounds. 
It is impossible to photoevaporate water (or any other 
high molecular weight material). Even if heated to 10 4 K, 
water cannot escape from the surface of the observed 
planets, because its sound speed would be ~ 2 km/sec, 
too low compared to the escape velocity^] Since the 
mass loss rate is determined by the smaller ofthe energy- 
limited and radiation/recombination-limited rates, and 
since the latter is exponentially suppressed when the 



sound speed is sm aller than the escape speed (Murray- 

inhib- 



■y 



Clay etaL 20091, the mass loss will be severe. 
itcdF 1 "! This conclusion differs from that of I Valencia et al.l 
(2010) who adopt the energy-limited mass-loss rate for 



water and hence find significant water loss. 

Another scenario, one that may have happened on 
Venus, is to photodissociate water molecules into Hi 
and oxygen and then remove H 2 via either Jeans escape 
or hydrodynamical outflow driven by photoevaporation. 
This would require a reductive agent to absorb the free 
oxygen. If this reductive age nt is atomic iron (see, e.g. 
Elkins-Tanton k, Seager|2008[ ), the original planet would 
have to contain 70% iron and 30% water in mass to ab- 
sorb all free oxygen into Fe^O^. Such a planet, at a mass 
of 10M e , would have a den sity of ~ 6g/cm 3 (c .f. the 
theoretical curves in Fig. [lj |Fortney et al. 2007), much 
denser than the cold planets we observe. In addition, for 

10 This sound speed may increase by a factor of a few if water is 
photodissociated and if both hydrogen and oxygen are ionized (E. 
Lopez, private communication). 

11 This argu ment means that the dus ty evaporative clouds from 
KIC 12557548 | |Rappaport et al.||2012| is hard to explain unless 
the surface escape velocity trom the planet is very low (E. Chiang, 
private communication). 



the newly produced hydrogen to be photoevaporated, it 
cannot be well mixed with the other heavy molecules. 
So we consider this option ineffective in compactifying 
Kepler planets. 

Instead we propose that the observed correlations can 
be explained if all compact planets started out as rocky 
cores overlaid with small amounts of hydrogen. Planets 
experiencing the strongest radiation would have their hy- 
drogen completely depleted, exposing dense rocky cores. 
We quantitatively investigate this possibility here. 

As a useful rule of thumb, at T cq = 1000 K, a hydrogen 
envelope of 1% in mass will expand the size of its planet 
by 25 — 40%, depending on the core composition and 
massH Another useful rule of thumb is the timescale to 
erode 1% of the planetary mass. We recast eq. (l8|) as, 



l%Mp 
M 



> 5Gyrs 



0.1 AU 



2.5i? fl 



7.5M 







M„ 



(9) 



where we have scaled the variables using typical param- 
eters for compact planets. Since typical planets in our 
sample orbit with period ~ 10 days, they just straddle 
the range where photoevaporation can make an order- 
unity change to their radii. 

For the compact planets, mass fractions in the hydro- 
gen envelope run from to ~ 1% (grey s ymbols in Fi. 



6 see assumptions for the calculations in £5.2.1). Fig 



also shows that smaller compact planets m our sample 
could have experienced evaporation loss that ranges from 
1% to 20% in mass; while larger compact planets typi- 
cally suffer loss of order 1% in mass or lower. This then 
makes a self-consistent story. Larger planets are larger 
because they can hold on to their hydrogen envelope, 
while smaller ones have lost it all, mostly as a result of 
their close proximity to the stars. 

We can stretch the calculation a little further. Assum- 
ing the planetary cores are not fully rocky (as assumed 
above) but are instead made up of a half-water/half-rock 
mixture, we would then require hydrogen mass fractions 
of order 10~ 3 or less to explain the observed planets. 
Such a thin atmosphere have no chance of survival for 
most of the planets in our sample. The situation is worse 
if the cores are made up of pure water. We conclude 
that the cores of these planets are likely composed of 
rock (or even denser material). This statement is sub- 
stantially weakened, however, if the internal luminosity 
of the planet is much lower than that of Neptune. 

So we are able to resolve the structural degeneracy in 
compact planets, helped by the fact that typical planets 
in our sample have T cq sa 1000 K. Our proposal, that 
compact planets are made of rocky cores overlaid with 
hydrogen < 1% in mass, explains the rock-like density 
of hot planets (e.g., Kepler-lOb, CoRoT-7b, 55 Cnc e), 
explains the factor of ~ 2 difference in size between hot 
and cold compact planets in the Kepler catalogue, and 
explains the density contrast within 2:1 pairs and the 
lack of density contrast in other closer pairs. 

12 However, the expansion is much smaller if the hydrogen at- 
mosphere is heavily polluted by metals. Photoevaporation is also 
problematic in that case. 
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The hydrogen mass fractions we infer for many of the 
compact planets are of order 1%. This can be primordial 
in orig in. Alternatively, as suggested by Rogers et al. 
pOlI ) and |Elkins-Tanton & Seager] ( |20'08f 7 outgassing 
and the subsequent break-down ot water could produce 
a hydrogen envelope < 1% in mass. The former option 
implies a fast formation timescale, while the latter im- 
plies a significant amount of primordial water (> 10% in 
mass) when the planet was formed. 



5.4. Rock or Water: the Lopez et al. (2012) study 



Lopez et al. (2012) study the thermal evolution of 
Kepler- lib, employing the energy-limited formula for 
photoevaporation. Their stellar FUV flux dims from 
6 x 10 4 erg/s/cm 2 at 100 Myrs, to 37 erg/s/cm 2 at 8 
Gyrs. When they consider the case that Kepler-llb is 
a rocky sphere overlaid with hydrogen, our preferred so- 
lution here, they find that it initially had to have thirty 
times more hydrogen and a much bigger size (for an ef- 
ficiency e = 0.02; see their Fig. 2). They also find that 
if the efficiency is 0.1, the initial planet mass would have 
been so large as to destabilize the intricate multi-planet 
system (Fig. 3 in that paper). They therefore conclude 
that Kepler-llb, with a density of 3g/cm 3 , is most likely 
a water-world. 



However, a concern with the Lopez et al. (2012) study 
is that they always adopt the energy-limited formula for 
mass-loss. If the planet is initially 8i?® and if e = 0.02, 
we find that mass-loss transitions from energy-limited to 
radiation/recombination-limited whenever the FUV flux 
exceeds w 10 4 erg/s/cm 2 . This transition occurs at an 
FUV flux 40 erg/s/cm 2 if e = 0.1. So the energy-limited 
formula over-estimates the mass-loss at the critical early 
stage when the star was FUV bright. Taking this into 
account may make the initial state of Kepler-llb more 
proba ble under the no- wa ter hypothesis than was made 
((2012). 



Lopez et al. 



out m 

To resolve these issues, it would be desirable to con- 
duct an evolutionary study of mass-loss that includes 
radiation/recomination-limited mass loss rates. We leave 
such a study to future work. 

5.5. Compact vs. Mid-sized 

We have suggested, based on the sharp drop-off in 
planet number at R = 3i? ffi , that compact and mid-sized 
planets are physically distinct categories. However the 
two types of planets span a similar mass range and ap- 
pear to form a continuum in their density and radius 
(Fig. [5]). It is unclear what induces one, but not the 
other, to accrete a healthy gas envelope, when both types 
of planets live in similar environments. We note that 
when a mid-sized and a compact planet coexist in a sys- 
tem, the former is always at least a factor of 2 or more 
massive than the latter. 

6. SUMMARY 

Based on two samples of low-mass planets — 22 planet 
pairs characterized by the TTV method and a dozen 
planets characterized by the radial velocity technique — 
we reach the following conclusions: 

1. TTV phases for most of our near-resonant pairs 
lie close to zero. We conclude that the majority 



of these pairs are consistent with an eccentricity 
distribution that has a root-mean-squared value of 
e ~ 0.01. About a quarter of the pairs, on the 
other hand, have eccentricities as large as 0.1 — 
0.4. True planet masses for the low-eccentricity 
population can be recovered statistically, based on 
their inferred eccentricity distribution. 

2. Masses of planets are roughly proportional to their 
radii, such that the best-fit solution corresponds 
to a population that has a constant surface escape 
velocity of 20 km/ s. 

3. Mid-sized planets (i? > 37?®) in our sample in- 
variably have such low densities that they have to 
contain substantial H/He envelopes. Their current 
rates of photoevaporation suggest that their masses 
and radii may have been limited by this process 
and that masses in their cores should be at least 
comparable to those of their gaseous envelopes. 



4. Mid-sized planets, some 23% of the B12 sample, 
show up exclusively around stars more massive 
than 0.8M Q . Perhaps relatedly, planets around 
lower mass stars tend to have smaller sizes. We 
do not have an explanation for these observations. 

5. Densities of compact planets (R < 3i?©) fall be- 
tween that of pure rock and pure water, with some 
even less dense than water. Planets with higher 
equilibrium temperatures tend to be denser and 
smaller. Planets in our sample fortunately strad- 
dle the region where photoevaporation could signif- 
icantly erode a 1% hydrogen envelope. This allows 
us to break the degeneracy in internal composition 
and show that it is possible to explain the observed 
correlations if these planets are mostly likely rocky 
cores overlaid with a layer of hydrogen that is < 1% 
in mass. These planets are likely not water worlds. 

6. The data suggest that 3i?© is a dividing line be- 
tween 'super-earths' (largely solid) and 'hot Nep- 
tunes' (extensive gaseous envelopes). The atmo- 
spheres on mid-sized planets were most likely ac- 
creted, while those on compact planets may have 
come from accretion or outgassing. 

We foresee a number of directions to pursue in the fu- 
ture. First, more pairs are needed. We are currently re- 
stricted to pairs in the immediate vicinity of resonances, 
and ones that orbit with periods 5 — 12 days. Transit 
timing data of a longer span or a higher quality would 
allow one to probe many more planet pairs. It will be of 
interest to see whether the density- T cq correlation (Fig. 
[8| extends to larger ranges. More pairs will also pro- 
vide a more accurate eccentricity distribution. At the 
moment, uncertainties in mass determination arise both 
from measurement errors in TTV amplitudes, and from 
statistical uncertainties when converting nominal mass 
to genuine mass. 

Our conclusion that gas-rich planets appear to have 
suffered significant photoevaporation helps explain why 
there are so few of them close to the star. It also pre- 
dicts that more of these gas-rich planets will show up 
at longer periods. The extended Kepler mission should 
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be able to resolve the issue. It is also of interest to fol- 
low self-consistent thermal and photoevaporation evolu- 
tion for these planets, in order to determine their initial 
states. It remains puzzling how such low-mass planets 
could have accreted so much hydrogen, so close to the 
star. Lastly, the origin of hydrogen on compact planets 
is an interesting question to pursue. If it was accreted di- 
rectly from the disk, these planets had to form before the 
disks disperse. If it was outgassed from broken-down wa- 
ter, these planets had to have water-rich interiors when 
they formed. 

Most of our planets, gas-rich or compact, have masses 
< 2OM0, a decade lower in mass than the Jovian plan- 
ets at the same period range. This mass, incidentally, 
is roughly the gap-opening mass in a disk with a scale 
height H/R w 0.03. This may be a valuable clue for 



planet formation. 

Moreover, other processes not considered in this study 
may be important. We argue that an icy object cannot be 
converted into a rocky ball. But we have not considered 
ice removal by bombardment of planetesimals or proto- 
planets. 

And lastly, we need to understand the source of eccen- 
tricity in multiple low-mass planet systems. 

We thank Doug Lin, Eugene Chiang, James Owen, 
Jonathan Fortncy and Eric Lopez for interesting discus- 
sions. YW acknowledges the hospitality of KIAA-PKU, 
where part of this work was performed. This research 
is supported by NSERC and the province of Ontario, as 
well as NSF grant AST-1109776. This study would not 
have been possible without the amazing work by the Ke- 
pler team. 
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Fig. 11. — TTV data and periodogram for 6 newly identified pairs of Kepler candidates. For each pair, the three panels show the inner and 
outer planet's TTV (top and middle panels), and their periodograms (bottom). All planet candidates in this sample show clear sinusoidal 
variations at the desired super-period, marked by the red line in the periodogram. In the top two panels, vertical arrows are times when 
the longitude of conjunction points at the observer; green sinusoids are the best fit sinusoids at the super-period; and fitted amplitudes 
a nd phases (and 68% confidence limits) are as listed, where the phase is relative to the vertical arrows. The fitting procedure is described 
in |Paper I| 



